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A  method  is  presented  for  the  formulation  of  unknown  param¬ 
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dynamic  systems.  <  In  contrast  to' other  available  techniques,  this 
method  requires  only  as  many  measurement  sensors  withinthe  field 
as  there  are  unknown  boundary  conditions.  Results  are  presented  for 
simulated  data  from  an  example  of  heat  conduction  with  radiation 
boundary 'and  for  experimental  data  from  a  cantilever  beam  with  a 
nonlinear  moment  at  the  boundary.  The  method  may  be:  applied  to 
partial  differential  equations  which  are  linear,  one -dimensional, 
and  have  known  time  invariant  coefficients.  :  The  nonlinear  boundary 
conditions,  are  specified  up  td  a  set  of  unknown  constant  parameters 
which  appear  linearly  in  the  boundary  conditions. 
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Aba  tract 

A  method  la  presented  for  the  formula¬ 
tion  of  unknown  parameters  in  nonlinear 
boundary  conditions  in  distributed  parame¬ 
ter  dynamic  systems.  In  contrast  to  ocher 
available  techniques,  this  method  requires 
only  as  many  measurement  sensors  within  the 
field  as  there  are  unknown  boundary  con¬ 
ditions.  Results  are  presented  for  simu¬ 
lated  data  from  an  example  of  heat  con¬ 
duction  with  radiation  boundary  and  for 
experimental  data  from  a  cantilever  beam 
with  a  nonlinear  moment  at  the  boundary. 

The  method  may  be  applied  to  partial  differ* 
entlal  equations  which  ore  linear,  one- 
dlmcnslonal,  and  have  known  time  invariant 
coefficients.  The  nonlinear  boundary  con- 
editions  are  specified  up  to  a  set  of  un¬ 
known  constant  parameters  which  appear 
linearly  in  the  boundary  conditions. 

I  n  t  rode,  t  i  on 

Distributed  system  identification  must 
address  Itself  to  four  interrelated  prob¬ 
lems.  They  are: 

I.  Measurement  Restrictions 

Since  the  system  is  distributed  In  space, 
measurements  can  be  taken  anywhere  within 
the  field  and  ac  the  boundaries.  i  wc»er, 
for  a  method  to  be  practical,  successful 
identification  should  be  accomplished  with 
a  minimum  number  of  sensors.  Furthermore, 
measurement  locations  in  che  field  must  be 
selected  such  that  nodes  of  the  dominant 
modes  of  the  system  are  avoided^)1. 

II.  Simulation  Method 

Closely  tied  to  the  measurement  problem 
above,  la  the  selection  of  a  simulation 
method.  For  purposes  of  either  digital  or 
analog  simulation,  the  infinite  order 
partial  differential  equation  must  he  re¬ 
duced  to  a  finite  order,  or  equivalent 
lumped  parameter,  model.  The  simulation 
method  selected  should  not  require  an  ex¬ 
cessive  number  of  measurements  and  if  im¬ 
plemented  digitally  should  be  fast  enough 
to  be  practical. 

’Numbers  in  parentheses  designate  References 
at  end  of  paper. 


III.  Performance  Criterion 

In  order  to  accomplish  ths  identifica¬ 
tion,  some  criterion  must  be  selected  to 
measure  the  error  in  the  method.  For  a 
practical  method,  selection  of  this  error 
again  must  not  require  sn  excessive  number 
of  sensors.  In  particular,  en  error  de¬ 
fined  over  the  entire  spatial  domain  la 
Impractical . 

IV.  Optimization  Technique 

Finally  an  optimization  technique  must 
be  selected  to  minimize  a  performance 
criterion  in  some  stable  fashion  by  Identl- 
fylng  a  "best"  set  of  parameters. 

Systems  considered  In  this  paper  srs 
those  which  can  be  adequately  represented 
by  known  one-dimenalonal,  linear,  station¬ 
ary  field  dynamics  but  have  boundary  con¬ 
ditions  which  are  nonlinear  algebraic 
equations  in  the  field  energy  variables. 

Heat  conduction  in  a  continuum  with  a 
radiation  boundary  condition  la  an  example 
of  such  a  system.  The  method  presented  In 
this  paper  identifies  the  unknown  parametezs 
in  these  nonlinear  algebraic  relations  be¬ 
tween  the  energy  variables  at  the  bound¬ 
aries. 

Development  of  the  Identification 
Procedure 

The  method  presented  in  this  paper,  as 
shown  in  Figure  1,  takes  advantage  of 
several  properties  of  distributed  parame¬ 
ter  systems.  In  simulating  partial  differ¬ 
ential  equations  one  can  separately  simu¬ 
late  the  field  dynamics  and  the  boundary 
conditions  and  then  combine  them  to  obtain 
the  interactive  response^).  Furthermore, 
since  internal  spatial  data  appropriately 
lccated  are  adequate  along  with  the  field 
model  to  describe  response  throughout  the 
entire  field,  one  can  use  information  at  a 
sensor  located  within  the  field  to  deter¬ 
mine  the  energy  variable  response  at  the 
unknown  boundary.  Thus,  the  appropriate 
energy  variable  can  be  simulated  at  the 
unknown  boundary  using  only  as  many  appro¬ 
priately  located  sensors  as  there  are  un¬ 
known  boundary  con’itlona.  Having  recover¬ 
ed  the  boundary  energy  variable  by  such 
simulation,  the  boundary  forcing  than  can 
be  expressed  as  a  function  of  only  the 
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i  .xnn.  I 


unknown  boundary  parameters.  If  the  un¬ 
known  boundary  condition  la  appropriately 
modeled,  this  boundary  forcing  will  be  a 
linear  function  of  the  unknown  parameters 
even  though  nonlinear  In  the  field  vari¬ 
ables.  For  the  full  power  of  the  proposed 
method  to  be  utilized  a  linear  function  In 
the  unknown  parameters  Is  necessary. 


3u  3q 

3t  3x 


(1) 


The  symbols  are  defined  In  the  nomenclature 
section  of  this  paper.  Since  the  heat  con¬ 
duction  is  assumed  symmetrical  only  half  of 
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Figure  1  Block  Diagram  of  the  Basic  Identification  Method 


After  forming  the  boundary  forcing  at  the 
unknown  boundary  as  a  linear  function  of 
the  unknown  paraaetevs,  an  error  is  defined 
at  the  original  measurement  location  be - 
twaen  the  actual  measured  response  and  the 
simulated  response  using  the  unknewn  bound¬ 
ary  condition  as  forcing.  An  ISE  perform¬ 
ance  Index  la  then  minimized  yielding  a  set 
of  linear  algebraic  equations  which  are 
solved  digitally  for  the  unknown  parameters 
(3.4).  Since  a  digital  computer  la  used  s 
steep  descent  technique  Is  net  necessary  to 
solve  the  optimisation  problem. 

The  method  Is  best  explained  bj  example. 

Application  of  the  Identification  Method 

Example  1.  Heat  conduction  vlth  a  radia¬ 
tion  and  forced  convection  boundary  con¬ 
dition. 

The  problem  considered  in  this  first 
example  Is  one-dimensional  heat  conduction 
in  «  homogeneous  slab  with  boundary  radia- 
ti  n  and  forced  convection.  Constant 
dlft  sivity,  environment  temperature,  and 
the  if.'tial  uniform  constant  temperature 
are  kn  .vn.  The  field  dynamics  are 


the  slab  la  simulated  and  the  boundary  at 
the  center,  x  »  0,  is  equivalent  to  an 
insulated  boundary.  At  the  unknown  bound¬ 
ary  condition,  x  »  1,  heat  trsnafer  occurs 
by  combined  radiation  and  forced  convection. 
Thus,  the  boundary  conditions  are 


q (0 , t)  -  0 


(2) 


q(l,t)  -  M(u(l,t>  +  i)'  -  ec*J 
+  n((u<l,t)  +  1)  -  oc) 


(3) 


Data  were  available  la  the  literature  for 
this  problem  with  n  •  0,  5  •  1.0'3)  and  for 
rj  •  1.0,  8  *  l.o(^).  These  data  in  the 
literature  were  generated  at  x  »  1  by  an 
Iterative  numerical  scheae  given  the  fl  end 
8  values.  By  further  simulation,  valuaa  of 
u(.S.t)  vers  generated  aa  the  measurement 
data  for  the  Id.  cification  method.  Fig¬ 
ure  2  shown  the  necessary  steps  in  applying 
the  Identification  msthod  to  this  example. 


By  applying  the  Laplace  transform  to  the 
partial  differential  aquation,  an  ordinary 
differential  equation  la  obtained  In  terms 
of  u(x,s)  and  x  where  s  In  the  Laplace 
variable.  If  the  actual  seasurad  Internal 
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Figure  2  Boundary  Condition  Identification  Procedure 


data  at  x  ■  x,  is  u»(t),  then  the  spatial 
conditions  used  to  recover  u(l,t)  are 
Equation  (2)  and 

*  u„(t)  (4) 

Using  these  spatial  conditions,  a  trans¬ 
cendental  transfer  function,  Cx,  Is  obtain¬ 
ed  relating  u(l,t)  to  uB(t) 

u(l.s)  -  Cjtu.ls)}  (5) 

*°  simulate  this  Cx  transfer  function,  de¬ 
fined  In  the  Appendix,  truncated  Infinite 
product  expansions  were  used  as  proposed  by 
Goodeonv5).  To  evcld  noise  probleas  due  to 
lead  network  eiaulatton  Involved  In  Gx,  a 
tenth  order  least  squares  curve  fit  was 
used  to  smooth  the  Input  data  without  Intro¬ 
ducing  phase  shift.  This  trsnsfer  function 
**•  then  digitally  simulated  using  e  matrix 
exponential  routine  with  a  tenth  order  hold 
on  t^fvinr*it  sa  developed  by  Krouse  and 
Ward*6). 

Having  recoversd  boundary  temperature, 
ud»t)»  the  unknown  boundary  hast  flux  Is 

formed  as 

gd. t)  -  eSjCt)  +  ns2(t)  <6) 

whsre  S1(t)  -  {(u(4,t)  +1)'  -  ©c'j  (7) 

and  S2(t)  -  ( (u(l,t)  +1)  -  9c]  («) 

Hcv  the  response  at  the  measurement  location 
can  be  estimated  as  a  function  of  the  un¬ 
knowns  8  andn  using  q(t)  as  s  spatial  con¬ 
dition.  Th-t  spatial  conditions  for  this 
simulation  tre  Equation  (2)  and  q(l,t)  In 
Equation  (6).  This  ylalds 


«(**»•)  -  C2lq(l,e)]  (9) 

etnee  q( 1 , t)  Itself  la  not  known,  but  Sx(t) 
and  S2(t)  ere  known,  the  estimate  of  uC5,  t) 
la  expreesed  se  a  function  of  the  unknown 
parameters  B  and  n. 

u(.5.s)  -  BC2!Sx(e)]  ♦  nG2l82(«)l  (10) 

To  identify  6  end  n,  an  arror,  e(t).  Is 
deflnsd  ss 

C(t)  -  R3(t)  -  IBRilt)  +  <l*2(t>]  (11) 

where 

i3(e)  -  um(s),  end  ijls)  -  CjfS^s)] 

1  -  1,2 

In  a  manner  similar  to  Kubin(T)  end  Kohr(®) 
an  ISE  performance  Index  Is  minimised  to 
find  the  "best”  5  and  n. 


(t)  dt 


Minimization  of  the  performance  index  re¬ 
sults  in  a  eat  of  linear  slgabraic  equa¬ 
tions  n  -  .  T  _  - 

6  rU  ■*  1  rl2  *  413 


*»  *22  "  *23 


where 


fM» 

M  ’  *1 

;‘s 


(t)  Ij(t)  dt 


Since  a  digital  computer  wee  used  In  this 
method,  this  set  of  linear  algebraic  equa¬ 
tions  was  solved  directly  using  s  Gaussian 
elimination  routine. 
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Identified  S  -  0.974 
u(.5,t)  BRjtt) 


>.60  -  ^ 

Estimated  radiation 

)a70  X  contribution  to  u{.5,t)./0 

Figure  3  Simulated  R^(t)  compared  with  Internal  Input  Temperature,  u^ft) 

x^  «  0.500 


U.ing  data  at  x„  -  0.5CJ  fro*  Croable 
and  Vlskanta  for  radiation  alone  at  the 
boundary  (n  -  0,8  -  1.0)  the  elmulated  re- 
aulra  in  Figure  3  were  obtained  aud  8  - 
0.974  was  identified.  Thus  the  radiation 
coefficient  was  identified  to  within  3X  of 
the  correct  value.  Using  this  identified 
value  of  8  the  agreement  with  the  actual 
measured  data  is  quite  good. 


In  Figure  4  the  resul.a  are  shown  using 
dota  at  xm  -  0.S00  for  combined  radiation 
and  forced  convection  at  the  boundary  (8  “ 
1.0,  n  “  1.0).  In  this  case  the  sensitivity 
to  the  radiation  contribution  is  much  less 
than  the  convection  contribution.  Thus, 
the  radiation  coefficient,  8,  was  identi¬ 
fied  only  to  within  20X  of  the  correct 
value  but  the  convection  coefficient,  n. 
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Figure  4  Simulated  Rj(t)  and  Rj(t)  Compared  with  Actual  Values  and  Total 
Input  Temperatura,  u  ft) ,  at  x_  -  0.500 
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was  correct  to  within  12.  More  importantly, 
the  estimated  u(.5,t)  resulting  from  these 
identified  parameters  virtually  coincides 
with  the  measured  data.  It  was  also  ob¬ 
served  that  a  higher  order  model  for  Gi 
would  allow  better  identification  of  3  (with¬ 
in  5-102)  because  of  inclusion  of  less 
significant  eigenvalue  terms. 

Engineering  judgment  must  be  used  to 
select  the  order  of  the  simulation  model. 

A  more  complete  and  therefore  more  costly 
simulation  model  results  in  the  more  accu¬ 
rate  quantitative  parameter  identification. 
However,  a  much  simpler  field  model,  while 
resulting  in  less  accurate  parameter 
identification,  may  still  yield  equally 
accurate  state  estimation  results.  In 
actual  practice,  the  agreement  between 
measured  and  estimated  responses  would  gen¬ 
erally  be  the  final  measure  of  a  successful 
identification.  Since  the  field  dynamics 
are  known,  well  known  linear  system  tech¬ 
niques  such  as  frequency  response  can  be 
used  to  determine  which  eigenvalues  are 
significant  and  therefore  where  the  model 
should  be  truncated.  A  more  complete  con¬ 
sideration  of  field  model  selection  is 
found  in  ref erence ( 2) . 

Example  2.  Transverse  beam  vibration  with 
a  nonlinear  support. 


Figure  6  Test  Beam  Mounted  in  Clamped- 
Unknown  Boundary  Configuration 


inertia,  is 


2il  ii 

EI  3t 


0 


(16) 


Equation  (16)  assumes,  as  a  first  approxi¬ 
mation,  that  an  equivalent  viscous  damping 
term  is  a  valid  approximation  for  internal 
structural  damping.  Since  the  system  was 
started  from  rest  both  initial  conditions 


In  the  second  example,  a  long  slender 
beam,  sketched  in  Figure  5,  was  excited  in 
transverse  motion  by  an  electrodynamic 
shaker  cable  clamped  to  the  right  end  of 
the  beam  as  shown  in  Figure  6  and  Figure  7. 
The  boundary  condition  at  x  -  0  was  the  un¬ 
known  nonlinear  one  which  was  identified. 
Furthermore,  the  energy  variables  at  this 
boundary  could  not  be  measured  directly  due 
to  the  interference  of  the  knife  edges. 

The  input  motion  of  the  shaker  cable.  z(t), 
was  measured  with  a  relucta.-ce  distance 
detector  and  the  measurement  internal  to 
the  field  was  obtained  from  a  set  of  stra'n 
gages  mounted  to  provide  s.  temperature  com¬ 
pensated  measurement  of  bending  moment. 

A  partial  differential  equation  model  foi 
this  system,  neglecting  shear  and  rotatory 


Figure  7  Beam  Test  Stand 


Figure  5  Test  Beam  Mounting  Configuration 
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*r*  zmrc .  The  known  spat  it1,  conditions  fit 

x  •  1  srs 

y(l.t)  -  r (t)  (17) 

*0<i  yj:  (l.t)  •  0  (18) 

At  x  •  0,  the  knife  r jgee  restricted  verti¬ 
cal  motion  such  that 

y(0,t)  -  0  (19) 

The  problem  is  to  formti  it r  ami  identify 
the  fuurth  remaining  be. -  ary  condition. 

Since  the  section  of  the  bean  Inside  of 
the  boundary  at  %i  -  0  was  only  10  per  cent 
of  the  total  beam  length,  the  displacescut 
of  the  tip  of  the  baaa  via  '«nuued  to  be 
proportional  to  the  slope  at  x3  -  0, 
3y/3xi(0,t).  Because  of  the  eventual  con¬ 
tact  vtth  the  rubber  stops,  the  force  _x- 
erted  ‘,n  the  tip  of  the  be-ij*  would  be  sia- 


ilsr  to  s  hardening  spring  fores.  Further- 
aore,  ths  aoaent  at  xj  -  0,  BI  3ly/3x2Z  (0,t), 
would  be  proportional  to  the  force  on  the 
besa  tip.  Thus  a  cubic,  hardening  spring 
relationship  was  essuasd  for  this  unknown' 
boundary  condition  model: 

M(t)  “  ^(0,t)  »  aiS,(t)  +  a2S2(t)> 

3*’  .  .,!)(.)  <J0! 

Ui“*  s1<t)  -  §*<0>t>  «l> 

s2(t>  -  (|^(0,t))2  (22) 

S3(t)  -  (f*(0,t))s  (23) 

The  &2  ter*  allows  for  a  nonsyaastrical 
boundary  1  ondltion.  This  boundary  modal  la 
able  to  fit  a  dead  band  reaaonably  well. 

Having  completed  the  important  modeling 
step,  the  application  of  the  identification 
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Figure  8  Simulated  y(.14,t)  Given  a'e  for  M  •  4, 
Cj  "  4Cj  Identification  Model. 

C2  ■  4c^  Damping  Model 
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method  proceeds  in  a  fashion  similar  to  « 
Example  1.  ,  Laplace  transforms  are  applied 
to  the  partial  dlffere-  ial  equation  in 
{16)  above.  The  three  known  spatial  con¬ 
ditions  and  the  internal  measured  bending 
moment  data,  Ua(t),  at  x  «  xB  were  used  to 
simulate  the  boundary  slope  at  x  -  0 
yielding 

|^(0, i)  -  G1l*(s)]  +  G2[ija(s))  (24) 

where  '“)  denotes  Laplace  transformed  vari¬ 
ables  ...-'d  again  Gi  and  C2  are  defined  in 
the  Appendix. 

Having  recovered  the  boundary  slope 
using  Equation  (24),  the  unknown  boundary 
condition,  M(t),  is  formed  as  a  function  of 
the  unknown  parameters  using  Equat ion  (20) . 
Using  h(t),  and  the  three  known  boundary 
conditions,  an  estimate  of  the  bending 
moment  at  x  -  xm,  p(xa,t),  is  simulated 
f  ton: 

0<*«.*)  *  C3lz(a)l  +  c4[M(s)]  (25) 

Now  the  error  at  x  -  xa  is 

£(»)  •  0a(s)  -  (a^ISj^s)]  +  a2G4lS2(s)J 

1  +  t»3c4rs3(s)]  +  G3[ z (s) 1 )  <26> 

Linear  regression  analysis  using  an  ISE 
performance  indejc  again  yields  a  set  of 
linear  algebraic  equations  which  are  solved 
for  the  desired  a’s  as  in  Example  1. 

Some  difficulties  were  encountered  during 
the  simulation  due  to  the  equivalent  viscous 
damping  model  selected  in  Equation  (16). 


Thus,  the  damping  ratio  of  the  second  mods 
was  increased  by  a  factor  of  four  over  the 
first  mode  damping  ratio  after  comparison 
,wlth  experimental  data.  After  the  boundary 
parameters  were  identified  a  separate  closed 
loop  simulation  was  used  to  check  the  re¬ 
sults  against  the  actual  measured  data. 
These  results  are  presented  in  Figure  8  for 
the  identified  bdundary  condition  which  is 
plotted  in  Figure  9.  Although  ringing  -ntill 
exists  in  the  simulation  model  due  to  cu 
imprecise  damping  model,  it  is  significant 
to  note  that  the  identification  method  was 
nevertheless  able  to  successfully  match  the 
peak  amplitudes  after  the'  first  two  peaks. 
These  peak  amplitudes  are  Important  since 
they  aid  in  predicting  failure  of  a  part. 
Distortion  of  the'  fi-rst  two  peaks  wss 
caused  by  FM  tape  recording  of  the  measured 
data.  :  ' 

Summary  and  Conclualona 

.  In  summary,  to  implement  the  method  pre¬ 
sented  fn  this  paper,  the  uaer 

1.  Modela  the  linear  field  with  the  known 
boundary  conditions. 

2.  Formulates  the  unknown  nonlln-er  bound¬ 
ary  condition  as  an  algebraic  relation¬ 
ship  linear  In  the  unknown  parameters. 

3.  Solves,  ky  simulation,  for  the  boundary 
variable  necessary  as  Input  In  (2) 
using  Internal  measured  data  for  the 
missing  spatial  -condition  in  the  model 
in  (1). 


(0)  *  10, 


in./(in.*/in.*) 


Y-  2.S10 
“j  -  -5.550 
°3  -  218.61 


a*  lAW J  OlJ mAOS  A&htA 


4.  Generates  the  unknown  boundary  output 
variable  as  a  function  of  the  unknown 
parameters  using  the  boundary  oodel  In 
(2)  and  the  results  of  (3). 


5.  Goodson,  R.E.,  "Distributed  System  Si*' 
ulation  Using  Infinite  Product  Expan¬ 
sions,"  Simulation,  Vol.  15,  Ho.  6, 
December  1970,  pp.  255-263. 


5.  Solves,  by  simulation,  for  the  response 
at  the  original  measurement  location  as 
a  function  of  the  unknown  parameters 
using  the  results  of  (4)  for  the  miss¬ 
ing  spatial  condition  in  the  model  in 
(1). 

6.  Identifies  the  unknown  parameters  using 
regression  analysis  applied  to  an  error 
between  the  actual  measured  data  and 
the  estimate  in  (5). 

Thie  parameter  identification  method  was 
successfully  applied  ro  a  heat  conduction 
example  using  numerically  generated  data 
from  the  literature  and  to  a  beam  equation 
example  using  experimentally  measured  data. 
In  both  cases  the  estimated  response  using 
the  identified  parameters  was  in  good  agree¬ 
ment  with  the  actual  measured  data.  Tn  the 
heat  conduction  example,  where  *  comparison 
with  the  correct  parameter  values  was  poss¬ 
ible,  the  correct  parameters  were  identi¬ 
fied  to  within  3  per  cent  In  one  case  and 
in  the  other  case  to  within  1  per  cent  on 
one  parameter  and  due  to  low  sensitivity 
10  to  20  per  cent  on  the  other  parameter. 

A  particularly  Important  result  la  the 
indication  that  a  single  measurement  not 
at  the  boundary  includes  enough  Informa¬ 
tion  for  identification  of  parameters  in 
the  boundary  given  some  knowledge  of  the 
field  dynamics. 
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APPENDIX 


Transfer  Functions  and  Inflnlea 
Product  Approximations 


Example  1  Heat  Conduction 
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For  xm  -  0.50C  and  a  fourth  order  mode) 
(N  -  4), 
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For  xtt  -  0.500  and  a  fifth  order  model, 
(N  *  4) , 
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Example  2.  Beam  Equation 

ylsinh  y  +  sin  y)(ainh  y  x  +  aln  y  x  )) 

Gi - r - 5 - - 

Denom. 
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For  x_  -  O.Ih  and  a  fourth  order  model, 
infinite  products  yleld(2) 
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c  .  Y3  (cos  Y  alnh  yx  +  cosh  y  sin  yx) 
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For  x  »  0.14  and  a  fourth  order  aodel, 
infinite  products  yields 
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Nomenclature 


“  internal  damping  term 
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transfer  functions  defined 
In  Appendix 
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boundary  parameters  to  be 
identified  in  beam  example 
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